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ABSTRACT 


We introduce two different formulae of mechanical energy of harmonic 
oscillators driven by an external force in this work. In the first formula, the energy 
is composed of two terms, i.e., the kinetic energy and quadratic-potential-energy 
terms. Another formula is that it is just the same as the Hamiltonian itself. While 
the first kind of energy reflects the mechanical energy remained in the system at a 
certain instant of time, the second energy depends on the linear potential term as 
well as the terms quadratic in canonical variables. Time evolution of the two 
energies is analyzed for two specific cases of which driving force is different from 
each other. The second kind of energy can be negative, whereas the first kind is 
always positive. The second mechanical energy has generally been used in 
optimal control theory in relation with shortcuts to adiabaticity. 


Keywords: Mechanical energy, Harmonic oscillator, Hamiltonian, External force, 
Classical solutions 


1. INTRODUCTION 


Mechanical energy is a basic element in describing various dynamical systems 
including oscillatory ones. Transfer of energy for a system interacting with other 
ones or its environment has attracted much interest in the physics community for 
a long time (Mukherjee and Barbatti, 2022; Chen and Muga, 2010; Choi, 2021; 
Hioe et al., 1978; Koochi and Goharimanesh, 2021). Mechanical energy is 
evaluated from the Hamiltonian of the considered system in most cases as we 
know. However, mechanical energy is not always the same as the Hamiltonian 
for complicated systems (Kobe, 1990; Kobe et al., 1990). For instance, in damped 
oscillatory systems, mechanical energy and the Hamiltonian are different from 
each other by a time-dependent exponential factor (Colegrave et al., 1989; Yeon et 
al., 1987; Choi, 2013). The role of the Hamiltonian in such a case is confined to the 
generator of the equations of motion of the system. The relation between energy 
and the Hamiltonian is delicate, while such a theme is not actively discussed 
previously. 

A notable consequence in this direction is that mechanical energy in driven 
oscillatory systems is defined in two ways in the literature. We will manage such 
two formulations of energy regarding their fundamental mathematical and 
physical aspects. Time evolution of the two energies will be analyzed for the cases 
where the type of the force is specifically chosen. From this, the differences and 
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similarities of the two energies will be clarified. We invoke readers the trend of such duplicate energy formulations on one hand, 
since most of researchers do not recognize it yet. 


2. MATERIALS AND METHODS 


We will look at the methodology of defining two different mechanical energies for oscillatory systems in this section. To analyze the 
physical characteristics of the two energies, it is not necessary to introduce a complicated system. We consider a forced oscillator of 
which classical equation of motion is simply given by 


X¥+a’x =F, (1) 


where q is the angular frequency and f(®) is an arbitrary force. The Hamiltonian that gives this equation of motion can be written 
as 
pi 


H =-—+-m[o’xX’ -2fOxl+e0, (2) 
2m 2 


where g ( is a time function which is real. The mechanics relevant to this Hamiltonian, including its solutions, is well known. The 
classical solution of position (momentum) is composed of a complementary function x, (p,) and a particular solution x, (p, ), 
such that 

XxX =X + Xp (3) 
Po =Po + Py. (4) 


The complementary functions for position and momentum are independent of the force and are given by 


xX, =X, sin@t+ 6), (5) 


D, =™M oX, cos (wt +6), (6) 


where x, is an amplitude and @ is an arbitrary phase. However, the particular solutions are different depending on the form of 
PQs 

Among the two classes of mechanical energy appeared in the literature for the above system, we first manage the one defined as 

— Bat + ym @ X2,, (7) 

where x,, and p,, are given in equations (3) and (4), respectively. This type of energy was claimed by Landau and Lifshits (1981), 
Dybiec and Gudowska-Nowak (2017), Romero-Bastida and Lopez (2020), Yeon et al. (2001), Um and Yeon (2002), etc. Yeon et al. 
(2001) and Um and Yeon (2002) treated mechanical energy of the forced damped harmonic oscillator, but, if the damping factor is 
removed from their treatment, their Hamiltonian and energy are reduced to equation (2) and equation (7) respectively. In 
particular, Romero-Bastida and Lopez (2020) used this concept of energy in analyzing the performance of an energy harvester 
driven by colored noise, which converts the power supplied by external noise into electrical energy. 

Another definition of mechanical energy is of the form 


e= Pai rer [o’x2,-2fOx,1+ 20, (8) 
2m 2 

which is naively the same as the Hamiltonian. This kind of energy for the same system was adopted by Kuzmin and Robnik (2007), 
Guéry-Odelin et al. (2019), Guéry-Odelin and Muga (2014), Ding et al. (2020), Zhang et al. (2021), etc. This energy has especially 
been applied to shortcuts to adiabaticity, whose main concern is to find fast routes of final results in an adiabatic change by 
optimally controlling system parameters (Guéry-Odelin et al., 2019). 

Let us call equations (7) and (8) as just the first energy and the second energy, respectively, for convenience from now on. The 
first energy is independent of g @) and the associated potential energy is measured from X = 0. However, the second energy is 
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different depending on g (6, and the reference point of the associated potential energy is different from that position. To see this in 
more detail, let us express equation (8) as 
2 


Pa 1 2 = 2 
ay tye ix,, -a@l-4@, (9) 


e@= 


where a@ = f@/a@* and b@ = mf’? @/@a’)- ¢@. We can think from this that potential energy in the second definition is 
measured from the reference position x = a (f) and — 6 @ is added in the energy as a compensation. 4(f) is used as a control 
function in shortcuts to adiabaticity, whereas 6 (¢) acts as a gauge function (see, for example, p. 37 of Ref. (Guéry-Odelin et al., 


2019). 


95» 1 a ar | 
(b) 20 | —— o,=0.5 —— o,=0.9 —— w=1.5 


(c) 20 —— 0,=0.5—— w,=0.9 —— wy=1.5 


0 20 40 60 80 100 120 


Figure 1 Mechanical energies with the choice of f() as equation (10) and g (f = 0. (a): time evolution of Xo with Xx, and Xyr where 
we have chosen q, = 0.9. (b, c): time evolution of energy where (b) is the one defined in equation (7) and (c) is the second energy e, 


. Panels (b) and (c) were drawn for several different values of g, as can be seen from the legend. The parameters that we used are 


o=1,x%,=1,§ =1,m =1, and @ = ¢ =0: These values are chosen to be dimensionless for convenience and this convention will 


also be used in the subsequent figures. 


3. RESULTS AND DISCUSSION 


Analyses with the choice of g (t) = 0 
We now see the two energies with the choice of g ({) = 0. Kuzmin and Robnik (2007) investigated how energy evolves when the 


oscillatory system is nonconservative with this choice of g @). Since the second energy is different depending on g (2), let us denote 


it as e, instead of e in this case. We consider two types of f(@ . At first, we put the force as a sinusoidal form such that 


FO =f sin@,t+ 9), (10) 


where h is the strength of the force and @ , isa driving angular frequency. While complementary functions are common regardless 
of the type of the chosen force, particular solutions are different depending on f() as mentioned earlier. The particular solutions 
associated with the force of equation (10) are given by 
4 , 
= sin@,t+ ¢), (11) 


p 
o -a; 
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f 
p= 71h cos (a,t+¢). (12) 
a — a; 


From Figure 1, we see that the evolutions of the two energies are somewhat different from each other. For the case of @, = 0.9, the 
scales of energies are very large compared to other two cases where w, = 0.5 and w, =1.5 . This outcome is due on account of the 
near resonance of the driving frequency @, with the natural frequency @. The scale of F is larger than e, in this case. Both 
energies fluctuate fairly regularly according to the variation of the oscillatory amplitude caused by the interference of the system 
with the external driving field. 

In particular, for the case of wm, = 0.9, the oscillations of energies are two types: the main oscillation with large amplitude (but 
with a small frequency) and the sub oscillation with small amplitude (but with a large frequency). The main oscillation accompanies 
the beat of x,,, which has taken place by the relative smallness of the difference between @ and q,. The sub oscillation has taken 


place on account of the direction of the exerting external force. That is, energy increases when the direction of the force is the same 
as the moving direction of the oscillator, whereas energy decreases when the direction of the force is opposite to the direction that 
the oscillator moves. Because such two relative directions occur alternately, the energies fluctuate with the sub-oscillatory angular 
frequency. 

The energies with @, =0.5 and @, =1.5 also evolve periodically but with complicated fluctuations. Such fluctuations also 


occurred depending on whether the direction of the force is the same or different compared to that of the oscillator motion. We can 
clearly confirm the difference in the evolution patterns of F and e for these two cases of O; from red and blue curves in Figure 


1(b, c), which are non-resonant with @. 


290 -—+ 1 1 1 a L 1 


(b) be —— a=0.1 == = a= a=0.3 


: \ ! . 
—— a=0.1—— a=0.2 —— a=0.3 


Figure 2 Mechanical energies with the choice of f(@) as equation (13) and g @ =0. (a): time evolution of y,, with x, and Xys where 


we have chosen a = 0.2. (b, c): time evolution of the two energies (b) £ and (c) e,. Panels (b) and (c) were drawn for several 


different values of @. The chosen scales of parameters are w =1, x, =1, § =1,m =1,and 6=0. 


We now take the force, as another example, in the form 
F® = fe, (13) 


where q is areal constant. The particular solutions in this case are given by 
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h -at 
=e, (14) 
p, = beat (15) 
ora 


whereas the complementary functions are the same as equations (5) and (6). We see from Figure 2 that the evolution of the two 
energies is different from each other in this case, especially at some initial instant of time. F oscillates due to the effect of the force, 
where amplitude of such an oscillation decreases with time. This means that the oscillator gives and takes energy with the external 
field that induces the force. Such an exchange of energy is gradually attenuated by the decrease of the force as time goes by. 
Oscillation of mechanical energy also appears when the oscillatory system is driven by a frictional force instead of the external force 
(Choi, 2021). Any force exerted on a system may induce a change of energy in the system in general. 

On the other hand, e, is slightly lower than zero around t=0, but it rapidly grows with the progress of time until it reaches a 
quasi-stable value that fluctuates a little. e, no longer varies soon after that. In summary, the fluctuations of both energies gradually 


disappear and reach a constant value as the force disappears exponentially. It is noticeable that e, can be a negative value as can be 


seen from Figure 2(c), whereas Fis always positive. 


For the case of g (t) = mf? (£)/(2w?) 
Let us see the case that g () = mf * (/ Qa”), which corresponds to b () = 0. The second mechanical energy in this case is reduced 


to 


2 


Pa , | 
6)= i +m a [x,,-a@F. (16) 


The subscript II in the symbol of this energy is inserted in order to distinguish it with the energy in the previous subsection. We can 
confirm that this energy is always positive, whereas the energy e, can be negative as we have seen previously. The energy in 


equation (16) was used by some authors in order to formulate a minimal-time transportation theory for a particle confined in a 


harmonic trap (Guéry-Odelin and Muga, 2014; Ding et al., 2020; Zhang et al., 2021). The time evolution of e,,is given in Figure 3 for 
both cases of f() represented in equations (10) and (13). From this figure, we see that e,,is larger than ¢,, but smaller than F in 


most instants of time. 


(a) 59 ——e,—e,—E 


1 L 1 L 1 1 1 1 1 1 1 


0 20 40 60 80 100 120 


(b) *° | —— 6 ——e,—— E 


Figure 3 Time evolution of e,, given in equation (16) and its comparison with e, and F for the case where f() is chosen as equation 


(10) for (a) and as equation (13) for (b). Panel (a) was drawn with the choice of w, =0.9, a=1, x, =1, 4 =1,m =1,and 9=¢=0. 


(b) was drawn with the choice of a =0.2, @=1, x, =1, § =1,m =1,and 9=0. 
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4. CONCLUSIONS 


We analysed two different formulae of mechanical energy for forced oscillatory systems. The time evolutions of the two energies 
were compared with each other under specific choices of the external force. This comparison was carried out for the first time as far 
as we know. The first energy is basic where it is composed of the kinetic and the potential energies that are given in terms of p’ 
and x’, respectively. Because square of canonical variables is positive, this energy is always positive. 

The second mechanical energy for a forced simple harmonic oscillator is a sort of mimic the Hamiltonian. Because the energy 
associated with the driving force can be negative, e, treated in Sec. 3 can be a negative value. On the other hand, e,, managed in the 
same section is always positive. Although we have used a basic system in our analysis, this research may make researchers rethink 
the description of energy and its relation with the external force from a fundamental physical level. 


Authors contributions 
Both authors contributed equally. 


Ethical approval 
Not applicable. 


Informed consent 


Not applicable. 


Conflicts of interests 
The authors declare that there are no conflicts of interests. 


Funding 
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) 
(No.: NRF-2021R1F1A1062849). 


Data and materials availability 
All data associated with this study are present in the paper. 


KEPERENCES AND NOES 


Chen X, Muga JG. Transient energy excitation in shortcuts to 


2014; 90(6):063425. 


adiabaticity for the time-dependent harmonic oscillator. Phys 8. Guéry-Odelin D, Ruschhaupt A, Kiely A, Torrontegui E, 
Rev A 2010; 82(5):053403. Martinez-Garaot S, Muga JG. Shortcuts to adiabaticity: 
Choi JR. Analysis of novel oscillations of quantized Concepts, methods and applications. Rev Mod Phys 2019; 91 
mechanical energy in mass-accreting nano-oscillator systems. (4):045001. 

Axioms 2021; 10(3):153. 9. Hioe FT, MacMillen D, Montroll EW. Quantum theory of 
Choi JR. Analysis of quantum energy for Caldirola-Kanai anharmonic oscillators: Energy levels of a single and a pair of 
Hamiltonian systems in coherent states. Results Phys 2013; 3 coupled oscillators with quartic coupling. Phys Rep 1978; 
(1):115-121. 43(7):305-335. 

Colegrave RK, Croxson P, Kheyrabady E, Ramjit UAT. Energy 10. Kobe DH, Aguilera-Navarro VC, Kurcbart SM, Zimerman 
of a nonconservative system in quantum mechanics. Physica AH. Exact phases and probabilities for a spin-1/2 particle in a 
A 1989; 161(1):118-127. precessing magnetic field. Phys Rev A 1990; 42(7):3744-3750. 
Ding Y, Huang T-Y, Paul K, Hao M, Chen X. Smooth bang- 11. Kobe DH. Invariance of the generalised Berry phase under 
bang shortcuts to adiabaticity for atomic transport in a unitary transformations: Application to the time-dependent 
moving harmonic trap. Phys Rev A 2020; 101(6):063410. generalised harmonic oscillator. J Phys A Math Gen 1990; 23 
Dybiec B, Gudowska-Nowak E. Underdamped stochastic (19):4249-4268. 

harmonic oscillator driven by Lévy noise. Phys Rev E 2017;96 12. Koochi A, Goharimanesh M. Nonlinear oscillations of CNT 
(4):042118. nano-resonator based on nonlocal elasticity: The energy 
Guéry-Odelin D, Muga JG. Transport in a harmonic trap: balance method. Rep Mech Eng 2021; 2(1):41-50. 

Shortcuts to adiabaticity and robust protocols. Phys Rev A 13. Kuzmin AV, Robnik M. Energy evolution in the time- 


Discovery 59, e18d1023 (2023) 


6 of 7 


DISCOVERY | RESEARCH ARTICLE 


14. 


15. 


16. 


17. 


dependent harmonic oscillator with arbitrary external forcing. 
Rep Math Phys 2007; 60(1):69-84. 

Landau LD, Lifshits EM. Mechanics, 3rd Ed. (Reed Edu. Prof. 
Pub. Ltd., Oxford, 1981), Sec. 22: Forced Oscillations. 
Mukherjee S, Barbatti M. A Hessian-free method to prevent 
zero-point energy leakage in classical trajectories. J Chem 
Theory Comput 2022; 18(7):4109-4116. 

Romero-Bastida M, Lopez JM. Efficient harmonic oscillator 
chain energy harvester driven by colored noise. Sci Rep 2020; 
10:14306. 

Um C-I, Yeon KH. Quantum theory of the harmonic oscillator 


Discovery 59, e18d1023 (2023) 


18. 


19. 


20. 


in nonconservative systems. J Korean Phys Soc 2002; 41(5): 
594-616. 

Yeon K-H, Kim S-S, Moon Y-M, Hong S-K, Um C-I, George 
TF. The quantum under-, critical- and over-damped driven 
harmonic oscillators. J Phys A Math Gen 2001; 34(37):7719- 
7732. 

Yeon KH, Um CI, George TF. Coherent states for the damped 
harmonic oscillator. Phys Rev A 1987; 36(11):5287-5291. 

Zhang Q, Chen X, Guéry-Odelin D. Connection between 
inverse engineering and optimal control in shortcuts to 
adiabaticity. Entropy 2021; 23(1):84. 


7 of 7 


